Anisotropic paramagnetism of monoclinic Na^T^Oy single crystals 
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Abstract 

The anisotropic paramagnetism and specific heat in Nd2Ti2C>7 single crystals are investigated. 
Angular dependence of the magnetization and Weiss temperatures show the dominant role of the 
crystal field effect in the magnetization. By incorporating the results from the diluted samples, 
contributions to Weiss temperature from exchange interactions and crystal field interactions are 
isolated. The exchange interactions are found to be ferromagnetic, while the crystal field con- 
tributes a large negative part to the Weiss temperature, along all three crystallographic directions. 
The specific heat under magnetic field reveals a two-level Schottky ground state scheme, due to 
the Zeeman splitting of the ground state doublet, and the ^-factors are thus determined. These 
observations provide solid foundations for further investigations of Nd2Ti20y. 
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I. INTRODUCTION 



Rare earth titanates (R^T^Or) are an interesting family of materials showin 



electric 



2, rich ferro- 

and magnetic properties jjj, Q]. The common crystal structures 1| for these 
materials are cubic pyrochlore for small rare earth ions (Sm 3+ - Lu 3+ ) and monoclinic for 
large rare earth ions (La 3+ - Nd 3+ ). The magnetic properties of the former have been 
explored extensively during the recent years [3j]. The subtle interplay among the various 
interactions including crystal field, exchange and dipole-dipole, coupled with geometrical 
frustration specific to the pyrochlore lattice j^J, has led to exotic magnetic phases such 
Q 



as spin ice [6|-l8| and spin liquid [9 
external stimuli such as pressure 
cesses have also been identified 12 



. These phases are extremely sensitive to doping and 



10l| and magnetic fields 



ll|. Interesting dynamic pro- 



13(. On the other hand, the rare earth titanates with 



a monoclinic structure have received much less attention. Due to the removal of geometric 
frustration, these materials are believed to possess magnetically ordered phases at low tem- 
peratures. However, the competition between the crystal field and exchange interactions 
(plus dipole-dipole) with close energy scales may lead to interesting and unusual dynamic 
freezing processes towards the ordered states. In this regard, it is essential to retrieve the 
magnitude of the crystal field, exchange and the dipole-dipole interactions, in order to study 
the magnetic freezing processes. 

Here we focus on the monoclinic Nd2Ti 2 0y (NTO). Some features of NTO, such as the 
ferroelectricity the optical magnetoelectric effect jl4| and the photo-catalysis effect 15], 
receive current interests in these materials. In this work, we study the dc susceptibility 
and specific heat of NTO single crystals. The anisotropic paramagnetism is investigated 
in detail. We found that crystal field interactions dominate the anisotropy, and exchange 
interactions J ex are estimated by combining the results of the diluted samples. The ground 
state of Nd 3+ spins splits into five Kramers doublets in the crystal field associated with the 
C\ site in P2i(4) space group. The magnetic field splits the ground state doublet into two 
singlets, as is evidenced by the specific heat at different fields, showing a two-level Schottky 
scheme with the energy gap proportional to the field. These studies lay solid ground for 
further investigation of this system under doping and external perturbations. 
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II. EXPERIMENT 



Single crystals of Nc^T^Oy and Nc^-zLaxT^Oy we re grown by the floating zone method 
in flowing oxygen. Seed rods for the single crystal growth were polycrystalline Nd2Ti 2 07 and 
Nd2_ a; La :r Ti 2 07 sintered by the standard solid state reaction in oxygen atmosphere, with 
starting materials of stoichiometric mixture of Nd2C>3, La2C>3 and Ti0 2 powder (Alfa Aesar). 
Oxygen atmosphere was used to avoid oxygen deficiency. The lattice constants were obtained 
by fitting the peak positions of the x-Ray diffraction pattern on polycrystals obtained by 
crushing Nd2Ti 2 C>7 single crystals. The single crystals were cleaved along a direction, and 
the in-plane b and c directions were determined by the Laue diffraction. Dc susceptibility 
were measured with a Quantum Design superconducting interference device magnetometer 
(SQUID MPMS-5) with fields of up to 5 T. Dc susceptibility at higher fields were mea- 
sured by the vibrating sample magnetometer option in a Quantum Design physical property 
measurement system (PPMS). The in-plane angular dependence of the magnetization were 
measured using a commercial sample rotator in MPMS, the diamagnetic background of the 
sample holder was carefully subtracted. For the susceptibility measurement along different 
crystalline directions, the alignment of the crystal along the direction of the magnetic field 
was maintained within the error of 5°. Specific heat was measured in a PPMS using the 
thermal relaxation technique. 

The as-grown NTO single crystals are of high quality and large size benefiting from the 
floating zone method. The samples are monoclinic in space group P2j(4), the diffraction 
peaks of the polycrystals obtained by crushing NTO single crystals are well indexed as shown 
in Fig. QJa), and resultant lattice parameters are a = 13.002 A, b = 5.468 A, c = 7.678 A, 
and a = 7 = 90°, /3 = 98.52°. Due to the layered structure, the crystal cleaves naturally 
along the (100) direction, as is shown in the XRD pattern in Fig. [1] (b). In Fig. [1] (c), 
the full width at half maximum (FWHM) of 0.09° of the rocking curve further indicates 
good crystallinity of the sample. The single crystals were cleaved and shaped into disks 
with typical dimensions of diameter ~ 1.5 mm and thickness ~ 0.2 mm. The disk shape 
mitigated the effect of demagnetization factor in the in-plane anisotropy of mag netization. 
However, we still corrected the demagnetization using the calculated value 16] N = 0.75 
for H || a and 0.12 for H || be respectively. 
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III. RESULTS AND DISCUSSIONS 



A. Anisotropic Paramagnetism 

In Fig. [2] (al)-(cl), we show x{T) at 1000 Oe with H || c and H || a for a single crystal 
and a polycrystal sample, respectively. The zero-field-cooled and field-cooled x{T) overlap in 
the entire temperature range in study (2 - 300 K), showing typical paramagnetic behavior. 
However, x{T) is clearly anisotropic when comparing the behavior of the single crystal with 
H || c and H || a, as well as that of a polycrystal. Above 150 K, the system follows the 
Curie- Weiss law with the effective moment of 3.44 - 3.49 yU#, close to the theoretical value 



of 3.62 n B (an experimental value [17j is 3.5 The resultant Weiss temperature 9w 

varies from 10.4 K for H || c to —128.9 K for H || a. $w for the polycrystal, which can be 
taken as the average effect from all the crystalline directions, is —42.1 K. The negative Weiss 
temperature of the polycrystal indicates an effective antiferromagnetic correlation between 
Nd 3+ spins. By comparing Fig. |2] (al) with (bl), we see that the in-plane magnetization is 
more than three times larger than the H || a case, this suggests that the Nd 3+ spins prefer to 
align in the be plane. Fig. [2](a2)-(c2) show the isothermal magnetization curves for H along 
different crystalline directions. As is expected for a paramagnetic system, M(H) shows no 
hysteresis, is linear at high temperatures, and shows convex curvatures at low temperatures 
due to the spin polarization process. The M(H) curves measured with the field applied 
along different directions again show strong anisotropy. 

Moreover, the in-plane M(H) curves themselves show anisotropy. Fig. [3] (a)-(d) show 
the angular dependence of the in-plane magnetization. A two-fold symmetry is found at all 
fields and temperatures. The normalized in-plane anisotropy M(6)/M m id, where M m id is the 



mean value of M(6), i.e., ^M(^), is shown in Fig. [3] (e). The magnetization along the c 
axis is roughly two times larger than that along the b axis. As we already noted that the 
Nd 3+ spins prefer to align in the be plane, the c axis is thus determined to be the easy axis 
of Nd 3+ spins. 

We also show l/x(T) with fields along different in-plane directions in Fig. [3](f). l/x(T) 
at different in-plane angles are parallel to each other at T > 200 K, as expected for param- 
agnetism with the same effective moments but different Weiss temperatures. The angular 
dependence of Weiss temperatures 9w is shown in Fig. [3] (g). 
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Magnetic anisotropy, in particular magnetization anisotropy in paramagnets were previ- 



ously reported in rare earth compounds 



18 



19] . The primary origin of this anisotropy is the 



crystal electric field. When an ion is placed within a crystal lattice, the degeneracy of the 
Hund's rule ground state multiplet can be lifted by the electric field from the surrounding 
atoms. In a rare earth compound such as NTO, the cystal field can overcome the external 
field and pull the total moment resulting from the spin-orbit coupling into a certain direc- 
tion. The anisotropic exchange coupling also contributes to the magnetic anisotropy. To 
separate the contribution of crystal field from exchange to the magnetic anisotropy, a series 
of diluted samples Nd2_; C La. r Ti207 were studied. Nonmagnetic La 3+ has very close ionic 
radius to that of Nd 3+ , it dilutes the Nd 3+ spins and keeps the lattice almost unchanged. 
Thus, in a good approximation, the crystal field does not change with dilution while the 
exchange interaction decreases systematically with increasing La 3+ content. Because the 
Weiss temperature of the system reaches as low as —150 K for Ndi.i 5 La .85Ti2O 7 , and even 
lower for samples with lower Nd 3+ content, which makes it difficult to realize a reliable 
fitting using the Curie- Weiss law with the M{T) data in the temperature range of 2 - 300 
K, only the data of Ndi.i5Lao.8sTi207 is presented here. The Nd 3+ content is determined 
indirectly by calculating the effective moment of the powder sample. By writing down the 
magnetic Hamiltonian, it is possible to extract both exchange and crystal field contributions 
to the angular dependent magnetization by assuming a constant crystal field contribution 
for the La 3+ diluted sample. Following the approach in Ref . [20] , the magnetic hamiltonian 
consists of two parts, the exchange energy 3f ex and the crystal field energy M'cf (the dipo- 



lar interaction can be assimilated in the exchange energy 20|)- The exchange interaction is 
anisotropic and can be expressed as 

where a, b, c denote the three crystalline directions, and J® x , J\ x and J£ x are the exchange 
constants for Nd 3+ spins along a, b and c directions, respectively. The crystal field part, for 
the low-symmetry C\ site of the Nd 3+ , is known to be 

^CF = Ep=2,4,6 EL-p B p°p ( 2 ) 

where Bp are crystal field parameters, and Op are Stevens equivalent operators. It is known 
that the high temperature susceptibility of a system with only one spin species can be 



expressed as [20] 



1 _ l/rp , 2i(J+l) , (2J-l)(2J+3) I?.V.ll \ 

X C' +L 3Affc lOfc _f/ 2 > /o\ 

+0(1). 

where C is the Curie constant, c is the percentage of the magnetic ions on the spin lattice, 
in Nc^-zLazT^Cv case, c = (2 — x)/2. N and \t are the exchange tensor and the crystal 
field tensor, and can be expressed by the parameters in eq. [I] and eq. [2] respectively. When 
the three axes of the crystal are of at least second-order symmetry (symmetry no lower 
than L 2 ), eq. |3]can be expressed with the parameters in eq. [1] and [2] along each crystalline 
direction as 



1 . - Xfrp J(J+1) jg , (2J-l)(2./+3) v / a\ 

Xa C \J- L 3 ex * 10fe ^ I 

here, a denotes the crystalline direction a, b or c, J" x are the exchange constants as already 
appeared in eq. [H and u aa is the a direction component of tensor u . In specific symmetries, 
u aa can be expressed explicitly with the crystal field parameters B^. The Weiss temperature 
is then 



_ - J(J+1) ia (2J-l)(2J+3) (r\ 
U W — C 3 J ex lOfc Uaa \ ' 

However, in the NTO case, the symmetry of Nd 3+ site is barely C\, an explicit expansion of 
V is thus difficult. Therefore only the exchange parameter J ex is extracted and the second 
term in Eq. [5] is left as a single parameter 8cf to be determined. 

The results of the Curie- Weiss fitting for diluted samples suggest that the contributions 
to the Weiss temperature from the exchange interactions 9° x = c "^ J 3 +1 ^ J" x and crystal field 

interactions 9^ F = — - — —tki: of opposite signs. The fitted Weiss temperature 6y/ } 

and the contributions from exchange interaction and crystal field 9 ex and 9cf m all three 
crystallographic orientations for NTO and Ndi.i5Lao.85Ti 2 07 are listed in Table I. With 
these Weiss temperatures, the exchange constants are determined to be J^ x = 7.1 ± 1.1 K, 
J\ x = 6.9 ± 1.1 K and J° x = 4.2 ± 0.7 K. It can be seen that the exchange constants are 
positive in all crystallographic orientations. This illustrates the dominant role of the crystal 
field contribution to the magnetic anisotropy. While the exchange interactions are always 
ferromagnetic, the apparent exchange deduced from the Weiss temperature in different di- 
rections can be either positive or negative. This is because the crystal field contributes a 
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large negative part to the Weiss temperature. As seen from Table I, 9cf is a couple of times 
larger in magnitude than 9 ex along a and h axes and comparable in magnitude along the c 
axis, making the Weiss temperature large and negative along a, b directions and small but 
positive along the c direction. 

It is also worth to note that eq. @]is independent of the magnetic field. Indeed, we see 
that in Fig. [3] (e), the normalized in-plane anisotropy M{9) / M m id is nearly constant for 
different fields at the same temperature. Only when spins are polarized significantly by the 
field, i. e., M(H) curve is not linear anymore, M(9)/M mi( i at different fields bifurcate, as is 
the case at 2 K. 

As we have discussed, the monoclinic NTO does not accommodate geometrical frustra- 
tion. However, with the absence of a magnetic transition down to at least 2 K, and the 
average Weiss temperature 9w ~ —42.1 K, the frustration index as defined by / = \9w\/T c 



(Ref. |2l||) of this system reaches about 21, and is even higher along a and b directions. The 
large frustration index may indicate an unusual ground state and interesting spin freezing 
processes. These investigations are currently under way. 

B. Specific Heat 

The magnetic ground state of A h/2 Nd 3+ on the C\ site in the P2i(4) space group should 
split into five Kramers doublets due to the crystal field. Calculations on the Nd 3+ in the 
same symmetry in other compounds indeed revealed the five doublets scheme 22)], while a 
detailed calculation on NTO is still lacking. The specific heat measurement provides more 
information on the magnetic ground state. Fig. [4] (a)-(c) show the magnetic specific heat 
of NTO for H || a, H || be and that of Ndi.i5Lao.85Ti 2 07 for H || a at different fields. The 
magnetic specific heat is obtained by subtracting the lattice contribution using the specific 
heat data of a nonmagnetic La2Ti207 single crystal. The zero-field specific heat in principle 
detects any magnetic phase transitions, and the magnetic entropy associates directly with 
the ground state degeneracy. In Fig. H] (a), the zero- field magnetic specific heat shows 
no anomaly for a magnetic transition, which is consistent with the paramagnetic behavior 
above 2 K. There is a progressive upturn as T decreases, which signals the presence of a 
large portion of the magnetic entropy below 2 K. The magnetic entropy is therefore difficult 
to integrate due to the limited temperature range. However, this problem is bypassed by 
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using the specific heat in a magnetic field. As shown in Fig. @] (a), C mag {T) in a field shows 
a broad peak with the magnitude of around 3.6 J/mol- K at low temperatures, and the 
maximum moves to higher temperatures with its magnitude almost unchanged as the field 
increases. These indicate that the peaks are Schottky anomalies associated with the ground 
state doublet splitting in the presence of a magnetic field. The two-level Schottky specific 
heat is, 

Cschottky = R(j:) 2 (l+e A / T ) 2 ^ 

where R is the gas constant, and A is the energy gap between the two levels. 

Eq. describes the data nicely as shown by the solid lines in Fig. H] (a). It can be seen 
in Fig. H] (d) that the fitting parameter A scales linearly with magnetic field. This further 
proves that the gap between the two Schottky levels are caused by the zeeman splitting of 
the ground state doublet. Using A = gfisH, the g factor for NTO with H along a is found 
to be 1.99. The magnetic specific heat is also studied in Nd2Ti 2 07 with H || be and in 
Ndi.i5Lao.85Ti207 with H || a. For NTO with H || 6c, the g factor is 3.10. The anisotropy 
in g factor again points to the presence of strong crystal field, and these g values should be 
corroborated with future crystal field calculations. 



IV. CONCLUSION 



In conclusion, we have studied both single- and poly-crystalline Nd2Ti2C>7 with angle- 
dependent dc susceptibility and specific heat. The anisotropic paramagnetism is observed 
and is mainly attributed to the crystal field effect. The easy axis is determined to be the 
c-axis in the be plane. The contributions to Weiss temperature from exchange interactions 
and crystal field interactions are separated. The exchange interactions are ferromagnetic, 
while the crystal field contributes a large negative part to the Weiss temperature, along all 
three crystallographic directions. This leads to large negative 8\y values along a and b axes 
and small positive 6w along the c axis. The 4 / 9 / 2 ground state multiplet of Nd 3+ splits 
into five Kramers doublets. The ground state doublet further splits into two singlets in a 
magnetic field, and this two-level ground state scheme is revealed by the magnetic specific 
heat, g factors show anisotropy too, again due to the crystal field. These studies provide 
solid foundations for further investigations of this system, including the magnetic ground 
state and the spin freezing processes. 
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Unit: K 


pa 


pa 

Vex 


pa 


°w 


Vex 


e CF 


°w 


6°ex 


6 CF 


Nd 2 Ti 2 7 


-128.9 


58.3 


-187.5 


-62.1 


57.0 


-119.1 


10.4 


34.3 


-23.8 


Ndi.i5Lao.85Ti 2 7 


-153.8 


33.7 


-187.5 


-86.3 


32.8 


-119.1 


-4.1 


19.7 


-23.8 



TABLE I: The Weiss temperatures 9w for Nd 2 Ti207 and Ndi.i5Lao.85Ti207 along a, b and c 
directions. The fitting errors are estimated to be within 3 K. The contributions from the exchange 
interactions and crystal field interactions are isolated as 9 ex and 9qf using eq. [5j 
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FIG. 1: (color online) (a) The x-ray diffraction pattern of the polycrystals obtained by crushing 
Nd2Ti207 single crystals. Most peaks are indexed in the monoclinic structure in the space group 
P2i(4), the index of some small peaks are not marked explicitly, (b) The x-ray diffraction pattern 
of a Nd2Ti20"7 single crystal. Crystals cleave naturally along the (100) direction, (c) The rocking 
curve of the (400) peak with FWHM = 0.09° indicates good crystallinity. 
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FIG. 2: (color online) (al)-(cl): Zero-field-cooled and field-cooled susceptibility vs temperature 
at 1000 Oe for H || c, H || a of a Nd2Ti2C>7 single crystal and a polycrystal sample, respectively. 
Insets are fittings by the the Curie- Weiss law. (a2)-(c2): Isothermal magnetization curves at 2, 10, 
150 K for H || c, H || a and polycrystal respectively. 
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FIG. 3: (color online) (a)-(d): The angular-dependent in-plane magnetization M{6) of Nd2Ti20"7 
at 2, 10, 100, 300K in the field of 0.5, 2, 5 T. (e): The angular dependence of the normalized 
magnetization M/M m ^, where M m j^ is the mean value of M(6). (f) The temperature dependence 
of the inverse susceptibility of a Nd2Ti2C>7 single crystal with the in-plane field along 8 = 33°, 
57°, 75°, 93°, 120°, 147°, 165°, 183° and 207°. (g): The in-plane angular dependence of the Weiss 
temperature by fitting 1/%(T) above 200 K in (f). Dashed lines labeled with b and c denote the 
in-plane crystalline b and c directions as determined by the Laue diffraction. 
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FIG. 4: (color online) Magnetic specific heat of (a) A Nc^T^Oy single crystal with H \\ a; (b) A 
Nd2Ti207 single crystal with H \\ be (the field is not applied intentionally along a specific in-plane 
direction); (c) Ndi.isLao.ssT^Oy with H \\ a. Solid lines in (a)-(c) are two-level Schottky fittings; 
(d) The energy gap A obtained by the two-level Schottky fittings. 
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